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ABSTRACT 



Aims. We study the time variations of the cycles of 20 active stars based on decades-long photometric or spectroscopic observations. 
Methods. A method of time-frequency analysis, as discussed in a companion paper, is applied to the data. 

Results. Fifteen stars definitely show multiple cycles; the records of the rest are too short to verify a timescale for a second cycle. 
The cycles typically show systematic changes. For three stars, we found two cycles in each of them that are not harmonics, and which 
vary in parallel, indicating that a common physical mechanism arising from a dynamo construct. The positive relation between the 
rotational and cycle periods is confirmed for the inhomogeneous set of active stars. 
Conclusions. Stellar activity cycles are generally multiple and variable. 
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1. Introduction 

In the middle of the last century it was realised that certain ob- 
served features of late-type stars could be explained by magnetic 
phenomena (Kron 1947), similar to those that had been detected 
on the Sun. The cyclic pattern of solar activity had already been 
known for more than 120 years (Schwabe I1843I I by the time 
systematic research of magnetic activity of late-type stars com- 
menced. In 1966 long-term monitoring of the relative fluxes of 
cores of the Call lines - thought to indicate the strength and cov- 
erage of surface magnetism through the enhancement of chro- 
mospheric flux in the line cores - of solar type stars began (cf. 
Wilson 1968), to search for stellar cycles analogous to the solar 
case, and that monitoring program continued for more than three 
decades. Eleven years after systematic monitoring had begun 
(Wilson 1978) published landmark results on the first decadal 
survey of chromospheric variability and stellar magnetic cycles. 

While Wilson was early in his observational work, a contem- 
porary comment on the possibility of the existence of stellar cy- 
cles appeared in an editorial note of Detre ( 197 1 ) to a paper pre- 
senting photometric observations of BY Dra by P.F. Chugainov, 
as follows: The continuous observation of Chugainov' s stars 
would be extremely important because diagrams like that on the 
opposite page may reveal the existence of cycles similar to the 
solar cycle in these stars. 

Studying the late type eclipsing binary XY UMa, Geyer 
d 19781 1 suggested an activity cycle to explain the outside-eclipse 
variability of the system. Durney & Stenflo ( I1972D presented a 
first theoretical forecast, namely, with increasing rotation rate, 
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cycle periods should be decreasing. Baliunas et al. (11996^ sum- 
marized knowledge on the dynamo interpretation of stellar ac- 
tivity cycles prior to 1996. 

As databases relevant to decadal stellar activity continued to 
lengthen and contain more stars, study began on the variability of 
stellar activity cycles. Multiple cycles, analogous to the known 
solar multicyclic variability, were recovered by Baliunas et al. 
( |1996bl l from the records of the Ca II index. Using these Ca in- 
dex data supplemented with photometric results of active dwarfs 
and giants Saar & Brandenburg ( 1999 ) studied time variability of 
stellar cycles and multiple cycles on the evolutionaly timescale, 
draw the distribution of Prot/Pcyc in the function of the Rossby 
number and compared the results with theoretical models. 

At the beginning of the current century, several photomet- 
ric records became long enough to study photospheric magnetic 
cycles on a relatively large sample of active stars. Using conven- 
tional methods Olah et al. ( 120001 120021 ) derived activity cycles 
and multicycles for most objects in the present investigation, 
and Messina & Guinan ( 120021 ) for six young solar analogues. 
Radick et al. ( 1998) investigated the relation between the Ca 
index and contemporaneous photometric measurements on 35 
Sun-like stars. Lockwood et al. ( 12007) studied the magnetic cy- 
cle patterns seen in the photospheres and chromospheres of 32 
stars primarily on or near the lower main-sequence. In both those 
last two papers it was found that on a decadal scale, younger 
stars decrease in brightness when their chromosperic activity in- 
creases, while the older, less active stars increase in brightness 
when chromospheric activity increases, as is the case for the Sun. 
The existence and relation of photospheric and chromospheric 
cycles are well established. 
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On the basis of the variabiHty of the 11 -yr solar cycle, which 
fluctuates between approximately 9 and 14 years, and the as- 
sumption that stellar activity shares similar properties to solar 
activity, one may presume that stellar cycles should also show 
multidecadal variability. The first attempt to follow changing 
stellar activity cycles was made by Frick et al. (1 19971 1, who de- 
veloped and applied a modified wavelet technique suitable for 
data with gaps, and found a variable cycle of one of the targets 
in the present investigation, HD 100180, from its Ca II index 
record. Subsequent work of Frick et al. (2004') and Baliunas et 
al. ([2006 ) used double-wavelet analysis of the records of stars in 
the Wilson sample to study interdecadal activity variations. 

Recent efforts have focused on methods to predict solar ac- 
tivity based on different, earlier observations and on dynamo 
theory. For critical reviews of those methods see Cameron & 
Schussler ( l2007b . and Bushby & Tobias (2007). The knowledge 
of the cycle pattern of as many active stars as possible may yield 
improved insight on the solar activity seen in the context of stel- 
lar activity. To that end, we have developed a method suited 
for study of multi-decadal variability from gapped datasets. The 
method is presented in detail by Kollath & Olah (2008, here- 
after Paper I), where it is applied to solar data and used to re- 
cover a complicated pattern of multi-scale variability of the Sun 
in the last few hundred years. In the current paper we apply the 
method to a sample of active stars with photometric or spectro- 
scopic records that extend over several decades. 

2. Observations 

2. 1 . The sample 

The most important criterion for the selection of stars for anal- 
ysis was a long record with few interruptions. In most of the 
cases gaps in the records are seasonal interruptions, and we re- 
quire that the gap does not exceed two seasons. Limited gaps are 
crucial because the method of analysis requires equidistant data, 
and therefore interpolation during gaps in the data (i.e., between 
observing seasons) is necessary. 

The stellar sample consists of 21 objects, of which the 
Sun is studied in Paper I. Most of the stars have photomet- 
ric records, and V833 Tau has a century-long record if aug- 
mented with photographic photometry; six stars have Ca II 
index records. The photometrically observed active stars are 
listed in Table [T] together with their rotational (or orbital) pe- 
riods, spectral types, vsini, and inclination. The last column lists 
the extent of the records, in years. Single stars are AB Dor, 
LQ Hya, V410 Tau, and FK Com (about its possible binarity 
see Kjurkchieva & Marchev 120051 1, and the rest (V833 Tau, 
EI Eri, V71 1 Tau, UZ Lib, UX Ari, HU Vir, IL Hya, XX Tri, 
HK Lac and IM Peg) are synchronised, close binaries. In the 
case of binaries Table [1] gives the orbital period and for single 
stars the rotational period. In addition we study six objects from 
the Wilson sample (Wilson, 1978): HD 131156A, HD 131156B, 
HD 100180, HD 201091, HD 201092 and HD 95735. In our 
analysis, those targets are considered to be effectively single 
stars, though HD 131156A with HD 131156B and HD 201091 
with HD 201092 form wide, visual binaries. Their rotational pe- 
riods, spectral types and length of the Ca index records are listed 
in Table |2] 

2.2. Multi-decadal photometry 

The photometric records contain all published material for each 
star in the sample, in V-colour. Those stars are well studied. 



Table 1. The stellar sample. I. Photometric observations. 



star 


rot. per." 


sp. type 


Vi'i'ni 


i 


time-base 




(days) 




(kms-') 


(°) 


(years) 


AB Dor 


0.515 


KOV 


91 


60 


18 


LQ Hya 


1.601 


K2V 


27 


65 


25 


V833 Tau 


1.788 


K5V 


6.3 


20 


20 




photograph 


ic + photometric data 




109 


V410 Tau 


1.872 


K4 


74 


70 


34 


EI Eri 


1.947 


G5IV 


51 


46 


28 


FK Com 


2.400 


G4III 


155 


60 


28 


V711 Tau 


2.838 


G5IV/K1IV 


41 


40 


30 


UZ Lib 


4.768 


KOIII 


67 


50 


17 


UX Ari 


6.437 


G5V/K0IV 


39 


60 


23 


HU Vir 


10.388 


KlIV/III 


25 


65 


17 


ILHya 


12.905 


LOIII/IV 


26.5 


55 


20 


XX Tri 


23.969 


KOIII 


21 


60 


21 


HK Lac 


24.428 


KOIII 


20 


65 


50 


IMPeg 


24.649 


K2III 


26.5 


70 


29 



The stellar parameters are from Strassmeier ( I2002t . and the references 
therein. 

" in case of synchronized binary the orbital period is given 



Table 2. The stellar sample. II. Ca index measurements. 



star 


rot. per. 


sp. type 


time-base 




(days) 




(years) 


HD 131156A =f Boo A 


6.25 


G8V 


35 


HD 131156B =fBooB 


11.1 


K4Ve 


35 


HD 100180 = 88 Leo 


14.3 


GOV 


34 


HD 201092 = 61 Cyg B 


35.5 


K7V 


35 


HD 201091 = 61 Cyg A 


36.1 


K75 


35 


HD 95735 = GI 411 


54.7 


M2V 


34 



with numerous papers describing the records. Recent summaries 
for the stars are as follows: AB Dor: Jarvinen et al. (I2005I I. 
LQ Hya: Berdyugina et al. (l2002b and Kovari et al. ( |2004] |. 
V833 Tau: Olah et al. dSOOTT l. V4 I0 Tau : Strassmeier et al. 
(1997b), EI Eri: Strassmeier et al. ( |1997bl l, FK Com: Olah et 
al. (.2006), V711 Tau: Lanza et al. (I2006I I and Strassmeier et 
al. (11997b |l, UZ Lib: Olah et al. (|2002] |. UX Ari: Aar um-Ulvas 
& Henry (120031) . HU Vir : Strassmeier et al. ( 1997^ . ILHya: 
Strassmeier et al. ( |I997bll, XX Tri: Strassmeier et al. (! I997b| l, 
HK Lac: Olah et al. (fT997] l. IM Peg: Ribarik et al. (2003). For 
several stars the records were updated through 2007-2008 using 
new photometry from the Vienna APT (Strassmeier et al. |I997a| l. 

For further details, see the references of the cited papers. 
Periods of sparsely-sampled data leading to large gaps, some- 
times seen in the onset of records, were disregarded. 

2.3. Ca II index measurements 

The Ca II measurements were collected at the Mount Wilson 
Observatory, with two sets of equipment on two different tele- 
scopes, from 1966 to 1978 and after, for details see Vaughan et 
al. ( 1978 ). The relative Ca emission is the flux ratio of two 0.1 
run passbands centred on the cores of the H and K lines and two 
2 nm passband in the nearby continuum. 
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FAP 




0.005 0.010 0.015 0.020 0.025 0.030 



A(noise)/A(signal) 

Fig. 1. False Alarm Probability that a structure appears in STFT 
when Gaussian noise, with different amplitudes, is added to the 
record of LQ Hya. Five different standard deviations were used 
and the results are plotted from left to right for cr„ = 0.01, 0.02, 
0.03, 0.04, 0.05 mags. Details are in the text. 



3. Method 

We apply a time-frequency analysis called short-term Fourier 
transform (STFT) to the photometric and spectroscopic records 
to study time variability of the activity cycles that occur in the 
sample of twenty stars. As stated earlier, the analysis requires 
equidistant data, and therefore we had to interpolate through 
gaps among observations. To interpolate, we used a smoothing 
spline. Before calculating the spline interpolation, we removed 
the rotational signal and the next four strongest components in 
the periodogram near the rotation period; the removal of the next 
four strongest signals addresses removal rotation signals that are 
non-sinusoidal in shape or other activity changes on similar time 
scales from the record. Because rotational modulation is sub- 
ject to change owing to differential rotation, removing the ro- 
tational signal (and the additional signals of similar timescale) 
was done separately for subsets in the record, typically around 
200 days long, which were then averaged for the spline inter- 
polation. We argue that such a procedure is necessary, because 
unevenly sampled rotational modulation may alter the average 
seasonal light level, and consequently cause a false signal in the 
analysis for multidecadal variability. This step is especially im- 
portant for stars with high amplitude rotational modulation of 
large amplitude, arising from, e.g., high axial inclination to our 
line of sight. The FWHM of the Gaussian used in the time series 
is usually ~180 days. Note, however, that in the actual calcu- 
lation we use filtering in the Fourier space instead of temporal 
filtering. The method is explained in detail in Paper I, where 
tests of the effects of seasonal gaps, rotational modulation and 
observational errors are also found. 

The effect of active region growth and decay, as discussed 
by Messina & Guinan (2003 ) has little influence on our results 
because we removed the rotational frequencies plus four other 
signals of similar timescale, thereby significantly reducing the 
amplitudes of the modulation due to rotation, in most cases. 
Moreover, cycles below about 1 .5-2 years are not considered sig- 
nificant, as was discussed in Paper I, and below. 

To check the effect of added Gaussian noise on false fea- 
tures in the time-frequency diagram, we performed a Monte- 
Carlo simulation. Here we present the results of that test on the 
record of LQ Hya as an example; data of other stars in the sam- 
ple are of similar quality. First Gaussian noise was added to the 
original observational record. We performed the same prepro- 
cessing (averaging and spline-smoothing interpolation) as for the 



standard processing of the actual records. The STFT of the dif- 
ference between the noisy and the original smoothed data was 
calculated, and the statistics of the maximum amplitude (largest 
peak or ridge) was estimated from 10000 cases. Fig.[T]shows the 
False Alarm Probability FAP(z), e.g. the probability that a struc- 
ture with amplitude larger than z = A{noise)/A(obs) appears in 
STFT owing to noise {A{noise), and A{obs) is the maximum am- 
plitude in STFT; i.e., z is the inverse of the signal-to-noise ratio). 
The test was performed with different standard deviations (cr„) 
of the added Gaussian noise. On the figure from left to right the 
curves belonging to o-„ = 0.01, 0.02, 0.03, 0.04, 0.05 mag ai-e 
displayed. In case of Gaussian noise ~0. 01-0. 02 mag - which is 
poorer precision than typical photometric precision - the proba- 
bility of false signals with amplitudes over 0.01th of the highest 
signal is less than 10% and fast decreases toward larger astro- 
physical amplitudes. 



4. Results of the time series analysis from 
photometric data 

From analysis of the records we found changing and multiple 
cycles for most of the stars. Below we present a short descrip- 
tion of the cycle pattern for each object. The records and STFT 
results are plotted in Fig|2] 

In the time-frequency plots darker colours mean higher am- 
plitudes. The signals are modified by different amplification fac- 
tors that are employed to make the smaller amplitude signals 
better visible (see Paper I for more). The amplification factors 
are given electronically in Table [3] It is not straightforward to 
estimate the significance of the signals - moreover, no indepen- 
dent measure can be given. We consider a cycle or cycles signif- 
icant if at least two of these three features are true: a.) the signal 
runs throughout the period of the observations, b.) it has high 
amplitude or c.) when two or more components are changing 
parallel. Cycles shorter than about 1.5-2 years are considered to 
be insignificant: we showed in Paper I that from datasets having 
yearly gaps such cycles cannot be recovered safely. On the other 
hand, a long-term cycle is not significant if its length is commen- 
surable to the time span of the observations, i.e., when the data 
do not cover approximately two full cycles. 

AB Dor. The record of AB Dor is not well sampled, the 
number of observations per year is low and gaps longer than 
two years occur twice. The only clear signal is a cycle ~3.3-yr 
The other, weaker signal of ~2-yr appears in the beginning of 
the record, and is not considered reliable, and possibly owes to 
undersampling. Variability of high amplitude is reported on the 
time scale of ~20-yr (Innis et al. l2008) . but the record is not long 
enough to verify it from our analysis. 

LQ Hya. This star is a fast rotating, single star, one of the 
most interesting targets. The light variation is well-sampled by 
the measurements and all gaps are less than one year long. 
Previously, cycles between 11.4-11.6 years, 6.8-6.5 years and 
2.8-3.2 years were determined by Olah et al. (2 000 2002) from 
observations spanning a shorter period, and more recently, cy- 
cles of 13.8±2.8-yr, and its harmonic 6.9±0.8 and 3.7+0.3-yr 
were found from long term photometry by Kovari et al. (2004)). 
The present analysis recovers two short cycles ~2.5 and 3.6-yr, 
the latter being stronger in the middle of the dataset. Another 
cycle ~7-yr is present, and that period increases continuosly to 
12.4 years with the same amplitude, while a small amplitude ~7- 
yr signal remains. 

V833 Tau. The short (20-yr) photometric record of this bi- 
nary was supplemented by archive photographic data (Hartmann 
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Fig. 2. Results of the STFT analysis. In the upper panel for each segment (i.e., for each star) the observational data (grey, blue in the 
electronic version), the data prewhitened with the mean rotational periods and its five harmonics for each data subset (black, red in 
the electronic version) and the corresponding spline interpolations (solid grey (blue) and black (red) lines) are plotted. The lengths 
of the data subsets are given in the top line of each segment. The lower panels show the time-frequency distributions, where darker 
colours denote higher amplitudes. The amplitudes are also modified by the amplification factors (see on-line TableO. Top: AB Dor, 
LQ Hya, V410 Tau. Middle: V833 Tau photometric data, V833 Tau photographic -i-pho to me trie data, EI Eri. Bottom: FK Com, 
V7 1 1 Tau, UZ Lib. Details are in the text. 



et al. I1981I I, thereby extending the dataset to 109 years. We 
analysed separately the well sampled short and more sparse, 
long records; the long one includes also the photometric data. 
The short photometric record shows two modulations, with the 
weaker ~2.2-yr long and the stronger ~5.2-5.5-yr. Both are 



present throughout the record, which agrees with our previ- 
ous results (Olah et al. (2000 2002). The rotational modula- 
tion has a small amplitude (few hundredths of magnitude) re- 
sulting from the low inclination of the star (^ 20°). This low 
amplitude of rotation modulation is an important fact when we 
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Fig. 2. (cent.) Results of the STFT analysis. Top: UX Ari, HU Vir, IL Hya. Bottom: XX Tri, HK Lac, IM Peg. Details are in the 
text. 



study the century long dataset from the archive photographic 
measurements, even though only 1-2 points/year are available in 
Hartmann et al. ( 119811 ). However, those data represent the stel- 
lar brightness, within the usual error of the photographic mea- 
surements of ~0. 1 mag, because the rotational modulation itself 
has a lower amplitude than this value. Thus, a cycle amplitude 
above ~0. 1 mag is well documented by the photographic record, 
and is considered real. A long-term modulation is present dur- 
ing the 109-yr period of observations, with a large amplitude of 
~0.9 mag. A long cycle of ~27-30-yr is also recovered. 

V410 Tau. The only T Tauri-type star in the sample, 
V410 Tau shows rotational modulation with a large amplitude 
(often ~0.5-0.7 mag). However, its mean brightness has re- 
mained within ~Q.2 mag during the 34-yr length of the record. 
The analysis reveals a small-amplitude cycle of ~6.5-6.8-yr, 
which abruptly changes to 5.2-yrnear JD 2449500 and afterward 
slowly decreases in period. The long-term trend seems to follow 
this decrease. Our result supports that of Stelzer et al. (2003), 
who found ~5.4-yr cycle in yearly mean magnitudes after 1990 
(-2448500). While Steltzer et al. remarked that the cycle had 
been out of phase in the earlier observations, we suggest a dif- 
ferent, longer, and slowly variable cycle length in the beginning 
of the record. 

EI Eri. In the recently available record covering 28 years, 
two short cycles of ~2.9-3.1-yr and ~4.1-4.9-yr are present. 
They smoothly and in parallel increase and decrease; the longer 
one has a high amplitude between JD 2446000-50000. A long 



cycle of ~14-yr with variable amplitude is also seen. Two cy- 
cles of ~2.4 and ~16.2-yr had earlier been found by Olah et al. 
( 120001) . Subsequently, using a longer record, Olah et al. (120021 ) 
confirmed the ~2.4-yr cycle and revised the longer one to ~12.2- 
yr. The first determination of the long cycle was 11 + 1 years by 
Strassmeier et al. ( 1997b J ; all determinations of the long cy- 
cle point to an approximate decadal cycle. Rotational periods 
of EI Eri in different seasons has been investigated in detail by 
Washuettl et al. ( 120081 ). and no correlation is found between the 
seasonal periods (or multiperiods owing to differential rotation 
of 2-3 active regions at different latitudes for most seasons) and 
the cycle pattern. 

FK Com. One dominant cycle is present in the record; it 
changes smoothly between 4.5 and 6.1-yr. This is consistent with 
the earlier result by Olah et al. (120061 ). who found quasiperiods 
of about 5.2 and 5.8-yr in the longitudes of starspots, based on a 
record of 18 years, between JD 2446800-53200 (1987-2004). 

V711 Tau. The 30-yr record has the best observational cov- 
erage among our sample. The observations themselves show two 
waves, of ~18 and 9-yr, and they are clearly defined in the time- 
frequency diagram. Apart from those two periods, a cycle of 
~5.4-yr is present throughout most of the record. A short, weak 
cycle of ~3.3-yr is also apparent from JD 2445000 onwards. 
Several papers give estimates of the cycles of this well-observed 
system, and all are in accordance with the present results: Henry 
et al. ( 119951 ) suggests a cycle of ~5.5+0.3-yr and some evidence 
for a long one of ~16+l-yr from photometry, Vogt et al. ( 119991 ) 
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reported short cycles of ~3.0+0.2 and ~2.7+0.2-yr for the po- 
lar spot area and for the low latitude spots from Doppler im- 
ages. Lanza et al. ( I2006I I using a longer record of photometry 
found a cycle of ~3-5-yr with variable amplitude and a longer 
cycle, ~19.5+2.0-yr Berdyugina & Henry ( 2007 ) using an inver- 
sion technique on photometric data between 1975-2006 found 
~5. 3+0.1 and 15-16-yr cycles. Taking into account all these in- 
dependent results both from photometry and spectroscopy, the 
cycle pattern of V7 1 1 Tau is well documented. 

UZ Lib. The cycle length slowly varies from 4.3 to 3.1-yr. 
Earlier Olah et al. (2002) had found a cycle length of 4.8-yr. The 
long-term signal near 15-yr is comparable to the length of the 
dataset itself, and thus cannot be established as a cycle. 

UX Ari. Its record contains a significant cycle of ~4.6-yr 
that deacreases to ~3.6-yrby the end of the record. A long-term 
signal seems to decrease in parallel with the shorter cycle, from 
~25 to ~15 years; however, that longer period timescale is com- 
parable to the length of the record, and is thus not considered sig- 
nificant. On the other hand, the fact that the long cycle changes 
in parallel with the short one strengthens its reliability. We note 
that Aarum-Ulvas & Henry (2003) report an activity cycle of 
25 -yr on UX Ari. 

HU Vir. The dominant cycle length of this star is ~5.5-6-yr, 
in agreement with the values found by Olah et al. ( 120001 12002] |. 
The shorter-period signals are considered to be insignificant. 

IL Hya. A long gap (of approximately 4 years) occurs in the 
beginning of the record, which makes the early part of the record 
unacceptable for our analysis. Hence, we considerably shortened 
the record for study here; compared to the previous studies of 
Olah et al. (2000, 2002, 2007), we are unable the study of the 13- 
yr cycle found earlier. However, we confirm the previous result 
of Olah et al. (120071 ) on the changing nature of the short cycle 
(from 3.5 to 4.3-yr in that paper). The present result shows that 
a ~4.4-yr cycle persists throughout the 20-yr record, and a short 
cycle, increasing from ~2.7 to ~3 yr also appears. 

XX Tri. Apart from a long trend, only a weak ~3.8-yr cycle 
is present with significant amplitude only in the second half of 
the record. 

HK Lac. The first cycle lengths of this star were estimated 
by Olah et al. ( |2000] ) as -6.8 and ~13-yr Frohlich et al. ( |2006] ) 
extended the dataset with measurements from Sonneberg Sky- 
patrol plates, and not only confirmed the previous results but also 
found a third cycle with a length of 9.65-yr The origin of the 
third cycle lies in the changing cycle lengths of HK Lac, as was 
first communicated by Olah (120071 ) . The very beginning of the 
record was omitted from the present analysis becasue of scarce 
sampling. The Sonneberg results start well before the photomet- 
ric measurements, thereby lenghtening the record to 50 years; 
the photographic measurements partly overlap the photometry 
and fill some gaps. The present results, in agreement with all 
previous ones, shows a cycle varying slowly between ~5.4 and 
~5.9-yr, together with the cycle increasing from ~10.0 to ~13.3- 
yr. 

IM Peg. One cycle of ~9.0-yr is present, in accordance with 
the ~10. 1-yr cycle found by Olah et al. (120021) . A long-term sig- 
nal of 18.9 years is not considered significant, because it is ten 
years shorter than the record. 

5. Results of the time series analysis from Ca II 
index data 

Studying stellar activity cycles was one of the main aims of the 
Wilson project (e.g. Wilson 119781 ) and from the Ca II index 



Table 4. Rotational periods from Ca II index records. 



Star rotation periods in days 





(1) 


(2) 


(3) 


(4) 


(5) 


HD 131156A 


6.25 


6 


6.31 




8.2 


HD 131156B 


11.1 


11 


11.94 






HD 100180 


14.3 


14 




38 




HD 201092 


35.5 


35 


37.84 


36 


45,51 


HD 201091 


36.1 


38 


35.37 


36 


44, 47.5 


HD 95735 


54.7 


53 






28.5 



(1) present paper, (2) Baliunas et al. dl996ab . (3) Donahue et al. ( 119961 ), 
mean periods, (4) Frick et al. j2004| (5) Baliunas et al. ( 12006b 



records, fundamental results on decadal variability have been 
published. Ours is not the first attempt to find not only cycle but 
also rotational periods from that dataset. However, we wished to 
compare, using the method here, the results from photometry of 
stellar photospheres with those from Ca II H&K spectroscopy, 
which reflects the variability of chromospheres. 

5.1. Rotational periods 

The first thorough analysis of rotational modulation from the Ca 
II index data was carried out by Baliunas et al. ( |1996al ). That 
work is based on records of 112 stars, among them are all the 
six objects of the present investigation. All rotational periods 
found by Baliunas et al. ( |1996al ) agree well with the present re- 
sults, which is not always the case in later studies of other stars 
in the sample. Fig. [3] shows the period analysis of the six stel- 
lar records, all of which were prewhitened with the long-term 
variability (trends and cycle periods). The results are presented 
in Table |4] together with periods found in the literature. In all 
cases typical patterns arising from Fourier analysis of non-stable 
signals are seen: e.g., differential rotation of the stars produces 
significant power at frequencies around the main rotational fre- 
quency of the amplitude spectra. 

We split the records into three nearly equal parts in time and 
repeated the period search. The results are different among the 
three sets, primarily because of the effects of differential rota- 
tion, but short lived structures appearing and disappearing at dif- 
ferent positions may modify the derived periods as well. The 
values of periods detected in the subsets of the records are given 
in the electronic Table |5] Taken into account these results plus 
those from the literature we estimate uncertainties in the mean 
derived rotational periods of at least a few tenths of a day in 
case of shorter periods, and about 1.5-2 days for stars with rota- 
tional periods of a few weeks. Seasonal periods determined by 
Donahue et al. ( 119961 ) show a larger range of periods than seen 
in our analysis because of a multi-seasonal average in our anal- 
ysis of several seasons combined as one subset. Baliunas et al. 
(1996a) and Donahue et al. (1996) used datasets 10 years shorter, 
and Frick et al. (2004 ) and Baliunas et al. (2006 ) the same length 
records as we use in the current paper. In Fig. [3](lower panel) we 
show two examples, HD 131156A and HD 131156B in which 
the amplitude spectra of the entire record and part of it yield 
different periods. Note that in the presence of strong differential 
rotation but without knowing the rotational profile and the incli- 
nation of the star, it is nearly impossible to derive a precise mean 
rotational period. 
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HD 9573b prewhiteiied dataset 



HD lOOiaO prewhitened dataset 



HD 201091 prewhiteiied dataset 



HD 201092 prewhitened dataset 





HD 131156A prewhitened. dataset 



HD 131156A prewhitened dataset (detail) 




-o.a - 

40C&0 








Fig. 3. Rotational periods from Ca II index data. Upper panel, from left to right: HD 95735, HD 100180, HD 201091, HD 201092. 
Lower panel, from left to right: HD 1131156A and detail, HD 131156B and detail. Each segment consists of three panels: top: 
data prewhitened with decadal-scale variations, middle: power spectra with the rotational period marked, bottom: data folded on the 
derived period. 



5.2. Time-frequency results 

Below we present the new results of the cycles of each star indi- 
vidually, and compare them with the previous results. 

HD 131156A = ^ Boo A. Variability on two time scales are 
present. The shorter cycle is ~5.5-yr and exhibits a high ampli- 
tude in the beginning of the record; the cycle length decreases 
to ~3.9-yr at the record's end. A longer-period variability is also 
apparent, with a characteristic timescale of ~1 1 yr in the middle 
of the record. Baliunas et al. (,1995) marked the decadal variabil- 
ity of the record as "var," which means "significant variability 
without pronounced periodicity, on timescales longer, than 1-yr 
but much shorter than 25-yr..", and this is in accordance with 
our results for the record that is ten years longer; also, Baliunas 
et al. ( I2OO6I 1 find a long cycle of 13.2 years. 



HD 131156B = ^ Boo B. Only one long-term periodicity is 
found, and only in the beginning of the record, approximately 
~4.3-yr, although a secular change is also evident. Baliunas et 
al. ([1995| remark as "long", i.e., signifcant variability is found, 
which is longer than 25 years. 

HD 100180 = 88 Leo. We derive a more or less constant cy- 
cle of ~3.5-yr, which persists nearly through the entire record, 
but whose amplitude is strong in the first half of the record 
and decreasing thereafter Another variable cycle of ~13.7-yr 
appears in the beginning of the record; the period decreases to 
~8.6-yr by the end of the record, in agreement with the earlier 
results of Olah et al. (2007). Our results in the beginning of the 
dataset are similar to those found by Baliunas et al. (1995), who 
found two cycles, ~3.56 and ~12.9-yr in the record shorter by 
ten years. Analyzing the full record, Frick et al. (|2004| l found 
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HD 95735 = GJ 411 P (rot. ) = 54.7 days dt: 240 
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Time-frequency distribution (STFT): 
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Fig.4. The STFT of the Ca II sample, panels and colour codes are the same as in Fig. |2] Top: HD 131156A, HD 131156B, 
HD 100180, bottom: HD 201092, HD 201091, HD 95735. Details are in the text. 



periodicities of ~3.7 and ~10.1-yr The short cycle is consistent 
with our result, the ~10.1-yr value may be an averaged mean 
value of the of the long variable cycle we detect. 

HD 201092 = 61 Cyg B. The record for this star also ex- 
hibits two activity cycles: one is ~4.7-yr and persists thoughout 
the record; the other has a timescale of ~10-13 years. Both peri- 
odicities have variable amplitudes. The long cycle was estimated 
in previous studies as 1 1 .4, 11.1 and 1 1 .7-yr by Baliunas et al. 
(fT995j . Frick et al. (.2004J and Baliunas et al. C2006J . respec- 
tively. 

HD 201091 = 61 Cyg A. The high amplitude cycle seen in 
the record of this star has a mean length of 6. 7-yr, which slowly 
changes between 6.2 and 7.2-yr. A shorter, significant cycle is 
found in the first half of the record with a characteristic timescale 
of ~3.6-yr Probably because of the changing and double nature 
of the main cycles, previous determinations are slighly different 
but still consistent with our result: Baliunas et al. (1995) derived 
7.3-yr, Frick et al. ( 120041 ) found 7.12-yr and finally Baliunas et 
al. ( 120061 ) estimated 5.43-yr However, the most interesting com- 
parison to our results for HD 20 1 09 1 is that of Frick et al. ( 1 997 , 
their Figure 5, bottom left panel, and Figure 6, lower panel). 
They derived a changing cycle between ~6.7 and 7.9-yr, while 
the record at that time was about ten years shorter than now. 
Comparing Figure 5 of Frick et al. (1 19971 ) to Fig. |4] (bottom, 
middle), it is seen that the results are similar, though, because of 
our higher frequency resolution we could resolve two cycles in 
the beginning of the record. 



HD 95735 = GJ 411. The stronger cycle of this star is 
~3.9-yr, which is slightly shorter (3.4-yr) in the beginning of 
the record. A longer, 11 -yr cycle is also present with a smaller 
amplitude. Baliunas et al. (1995 ) classify the record as "var" 
(see the explanation at HD 131156A). The long cycle was es- 
timated as ~1 1.3-yrby Baliunas et al. C2006) . The rotational pe- 
riod found by Baliunas et al. (120061) is 28.5-d and by us 54.7-d; 
the shorter result possibly arises from the non-sinusoidal nature 
of the rotational modulation caused by two distinct activity cen- 
tres at spatially separated stellar longitudes. 



6. Discussion 

6.1. Cycle lengths 

The results of our analysis are summarised in Table|6] where the 
cycle lengths (yr, and in the electronic Table |2l in d) are listed 
in order of increasing value. In cases of changing cycle periods, 
arrows show the direction of the change between the extreme 
values (increasing, decreasing, or both). An 'L' denotes a multi- 
decadal variation, which may occur in addition to the detected 
cycle(s); however, a timescale cannot be given because of the 
insufficient length of the record to resolve the longer period. 

Most of the stars show multiple cycles. Those stars for 
which we list only one cycle usually show variability on longer 
timescales. The lengths of the records are too short to verify ad- 
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Table 6. Derived activity cycles. 



star Rotation period" Cycle lengths 

(days) (years) 



Ad Uor 


n c 1 c 
U.jl J 


3.3, L 
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1.601 


2.3, 3.0, /.U— 
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V OJJ iclLl 


i . / oo 




\TA 1 n Toil 


i.O /Z 


4.1<— J.Z, O.J— >D.o 


rii Jiri 


1 Q/n 

/ 


Z.y<-> J. 1 , 4. 1 <-*4.y, 14. U 




/inn 
Z.4UU 


/I 1 T 
4. j<->0. 1 , Lf 
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V / 1 1 lau 


Z.5 JO 


aa c: A qq- 17Q 

J.J, J. 4, o.o< — 1 / .y 


T TV T ;k 


4. /OS 


a 1 ,i A Q T 

j.l^ — 4. J, L, 


UX Ari 


6.437 


3.6«-4.6, L 


HU Vir 


10.388 


5.7, L 


ILHya 


12.905 


2.7^3.0, 4.4, L 


XXTri 


23.969 


3.8, L 


HK Lac 


24.428 


5.4^5.9, 10.0-^13.3, L 


IMPeg 


24.649 


9.0, L 


HD 131156A 


6.25 


3.9^5.5, 11.0 


HD 131156B 


11.05 


4.3, L 


HD 100180 


14.32 


3.5, 8.6<-13.7 


HD 201092 


35.53 


4.7, 11.5 


HD 201091 


36.06 


3.6, 6.2^7.2 


HD 95735 


54.7 


3.4^3.9, 11.0 



The cycle lengths, given in years, are listed in order of increasing value. 
For changing cycle periods, arrows show the direction of the change 
between the extreme values. 'L' denotes multi-decadal variations. 
" in the case of a synchronized binary, the orbital period is given 



ditional, longer cycles for AB Dor, HU Vir, XX Tri, IM Peg and 
HD 131156B (see Table©. 

The cycle lengths are not strictly constant or randomly vari- 
able, but show systematic changes during the period of obser- 
vations. We see two kinds of systematically changing cycles. 
The first is characterized by a cycle length that oscillates be- 
tween a lower and an upper value, and the second by a cycle 
length exhibiting a secular variation, i.e., an increase or a de- 
crease. However, these continuously increasing or decreasing 
cycle lengths may reverse over a longer observational period 
than currently available. 

In Paper I we studied the multi-decadal variability of the 
solar Schwabe and Gleissberg cycles during the last 250 years 
from Sunspot Number records. In Figure 2 (bottom panel) of 
Paper I we used the same method (STFT) as for the stars in 
the present paper. The results look similar - i.e., for the Sun, 
one cycle (Schwabe) varies between limits, while the longer one 
(Gleissberg) continually increases. However, the results for a 
longer time displayed on Figure 8 in Paper I indicated that the 
Gleissberg cycle of the Sun also varies between limits on a scale 
of ~50-200-yr, but to confirm the variability of the Gleissberg 
cycle required a longer, 500-yr record. By analogy from the anal- 
ysis of the longer solar record, the presence of a long-term trend 
may suggest an increasing or decreasing multidecadal cycle that 
is presently unresolved in the stellar records of short duration. 

Another feature of interest seen in our analysis of cycle 
lengths in stellar records is that of the dominant (highest am- 
plitude) cycle rapidly switching to another one. In EI Eri in the 
beginning of the dataset the long, ~14-yr cycle is dominant, and 
the short cycle of about 4.5-yr is only weakly present. Around 
JD 2447000 the amplitude of this short cycle increases consider- 
ably, and that high-amplitude cycle persists for more than 12-yr 
(~3 short cycles). Later the longer cycle again dominates. (The 



beginning is uncertain because of sparse data, but the end of the 
record is well sampled.) A similar pattern is seen in LQ Hya. 

Interesting feature is found in V410 Tau, EI Eri and UX Ari: 
two cycles in each of them which are not harmonics of each 
other, are changing parallel. This finding points towards the com- 
mon physical origin of these cycles, i.e., that both cycles are pro- 
duced by the same dynamo. 

6.2. Relations 

The cycle lengths of two stars were excluded from the fur- 
ther analysis. The first is the shortest period, ultrafast rotating 
(f „, = 0.459 d). Ml -2 dwarf star FY Dra, which might show 
a cycle of around 1 year from observations made every two 
months with gaps of only 4 months (see Vida 2007), but that 
cycle is uncertain. The other is the star with the longest rotation 
period in our sample, HD 95735 (P,o, - 54.7 d), a dwarf of spec- 
tral type M2, which is rotating -120 times slower than FY Dra, 
but has the same spectral type. Those two stars are the only M 
dwarfs which have reported activity cycles, and have the short- 
est (FY Dra) and longest (HD 95735) rotation periods, which are 
plotted in Fig. |5]and Fig. |6] but are excluded from the fits. Their 
internal structure may differ significantly from the other stars in 
the sample. 

Graphical representation of the results are given in Fig.|5]on a 
log - log scale of Pcyc/Pmt vs. 1 /Proi, as in the paper of Baliunas 
et al. (|1996bl Figure 1). The cycle period is expected to scale 
~ D' where t is the slope of the relation and D is the dynamo 
number. The slope for the results in Baliunas et al. ( 1996b ) was 

0. 74 (dotted line in Fig.|5l), and we get 0.81 + 0.05 for all the re- 
sults from the present paper and from Messina & Guinan ( 120021 1. 
and Frick et al. (120041 1. The fit to the shortest cycles determined 
by us gives a slope of 0.84 + 0.06. 

The log - log representation of the cycle period normalised 
by the rotation period gives a direct relation to the dynamo num- 
ber In comparison, the direct representation of the cycle vs. ro- 
tational periods in Fig.|6]has information about the tightness of 
the relation that is in Fig.|5]suppressed a little by the normaliza- 
tion and the scale. All symbols of Fig.|6]are the same as in Fig.|5] 
The Gleissberg timescale of the Sun is not plotted, because even 
its shortest value (50-yr, i.e., about 18000-d) would compress 
the figure along the y-axis. The slope of the relation for all the 
results of the present paper together with the values from the lit- 
erature is 39.7 + 10.6, whereas for the newly determined shortest 
cycles is 27.1 + 12.1, i.e., the relation is significant for all the data 
considered together, and marginal when only the shortest cycles 
are considered. 

However, an important warning should be recalled from 
Paper I: cycle periods of 1-2-yr recovered from records with 
approximately yearly gaps are highly uncertain if not insignif- 
icant. Hence, the lower end of the diagrams in Fig.|5]and Fig.|6] 

1. e., at the shorter rotational periods, the shortest cycles simply 
cannot be recovered with certainty, thereby distorting the fitted 
relations. 

Finally, we note that the stars studied in this paper and 
those from the literature that are plotted together in Fig. |5] and 
Fig. |6] represent a wide range of objects in the context of spec- 
tral types, characteristics of the observational data (photometry, 
spectroscopy, lengths of the records), and that the sample con- 
sists of both single stars and stars in binary systems. Therefore, 
while we find a general trend between the rotational periods and 
cycle timescales, a much tighter correlation cannot be expected. 
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-1 -0.5 
log(l/P,J 

Fig. 5. The quantity logiPcyd Proi) ^^-logil /P,-oi), as in Baliunas 
et al. ( |1996bl l. The plot includes the results of the present inves- 
tigation, and also results from the literature: Messina & Guinan 
( 1200 2' circles), and Frick et al. (2004, crosses). The results from 
the current photometric investigation are marked by filled dots, 
and filled triangles denote the results from the Ca II records 
(HD 95735 is excluded from the fits and the mean of its chang- 
ing cycles are plotted). The different cycles of the Sun are repre- 
sented by stars. A hexagon in the far right shows the position of 
the fast rotator EY Dra (excluded from the fits). The dotted line 
shows the original relation from Baliunas et al. ( |1996bl l, dashed 
line is the fit to all results, including those from the literature, 
and the solid line is the relation for the shortest cycle lengths 
determined by us. 
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Fig. 6. The rotational period-cycle period diagram, symbols are 
as in Fig. |5] The dashed line is the fit to all results including 
the ones from the literature, and solid line is the relation for the 
shortest cycle lengths determined by us. 



future only trends and tendencies can be suggested. Further de- 
tailed theory is necessary to model the complicated cycle struc- 
tures we observe in stars and in the Sun. Thorough modeling the 
observed cycle patterns may be for the present the only way to 
predict with success the future patterns of the Sun and the stars. 
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7. Summary 

We analysed continuous photometric and spectroscopic records 
on multi-decadal periods of 20 active, cool stars using Short- 
Term Fourier Transform and studied the patterns of their activity 
cycles. We found that 

- 15 stars show definite, multiple cycles, the others probably 
do, 

- most of the cycle lengths are found to change systematically, 

- in three cases two cycles are changing parallel, 

- and the relation between the rotational periods and cycle 
lengths is confirmed 

for the inhomogeneous set of active stars consisting of 
dwarfs and giants, single and binary systems. 

Moss et al. ( |2008l l depict a complicated picture of symmetry 
properties derived from models of stellar magnetic fields, while 
expecting that properties of observed stellar cycles may offer 
progress on the problem. The present paper gives twenty differ- 
ent examples (21 together with the Sun from Paper I) of multiple 
cycles that can arise from "dynamos operating in significantly 
supercritical regimes" . The observed modulations in amplitude 
and cycle lengths may originate from cyclic quadrupolar com- 
ponents superimposed on a dipole cycle, as suggested by Moss 
et al. ([2008). The predictability of the cycle patterns at the mo- 
ment is very limited: we are able to describe from observations 
with mathematical tools what happened in the past, but for the 
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Table 7. Derived activity cycles. 



star 


Rotational period" 


Cycle lengths 




(days) 


(days) 


AB Dor 


0.515 


1200, L 


LQ Hya 


1.601 


909, 13 16, 2571-»4545 


V833 Tau 


1.788 


794, 1898, 10000, L 


V410 Tau 


1.872 


1493<-1887, 2381^2500 


El Eri 


1.947 


1066«-^1136, 1506<-» 1786, 5100 


FK Com 


2.400 


1639<-^2222, L 


V711 Tau 


2.838 


1200, 1957, 3195<-6536 


UZ Lib 


4.768 


1130«-1560, L 


UX Ari 


6.437 


1299<-1667, L 


HU Vir 


10.388 


2100, L 


ILHya 


12.905 


993^1111, 1618, L 


XXTri 


23.969 


1400, L 


HK Lac 


24.428 


1957^2150, 3650^4854, L 


IMPeg 


24.649 


3298, L 


HD 131156A 


6.25 


1429^2000, 4000 


HD 131156B 


11.1 


1563, L 


HD 100180 


14.3 


1285, 3125<-5000 


HD 201092 


35.5 


1724, 4200 


HD 201091 


36.1 


1320, 2262<^2632 


HD 95735 


54.7 


1250<^1429, 4000 



The cycle lengths, given in days, are Usted in order of increasing value. 
For changing cycle periods, arrows show the direction of the change 
between the extreme values. 'L' denotes multi-decadal variations. 
" in the case of a synchronized binary, the orbital period is given 
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Table 3. Amplification factors for the different frequency intervals. 



star 


frequency (1/d) 


ampl. 


frequency (1/d) 


ampl. 


frequency (1/d) 


ampl. 


start-end J.D. 




interval 


factor 


interval 


factor 


interval 


factor 


-2400000 d 


AB Dor 


0.0 - 0.0004 


0.7 


0.00041 - 


4.0 






43856-50438 


LQ Hya 


0.0 - 0.0005 


1.0 


0.00051 - 


2.5 






45275-54450 


V833 Tau 1. 


0.0 - 0.0004 


1.0 


0.00041 - 0.001 


2.0 


0.0011 - 


5.0 


47098-54450 


V833 Tau 11. 


0.0 - 0.0004 


0.8 


0.00041 - 0.001 


5.0 


0.0011 - 


9.0 


14643-54450 


V410 Tau 


0.0- 


1.0 










41675-54011 


ElEri 


0.0 - 0.00035 


1.0 


0.00036 - 


3.5 






44129-54450 


FK Com 


0.0 - 0.0002 


1.0 


0.00021 - 0.0006 


2.5 


0.00061 - 


4.0 


43944-54285 


VVllTau 


0.0 - 0.0004 


1.0 


0.00041 - 


3.0 






42719-53703 


UZLib 


0.0 - 0.0004 


1.0 


0.00041 - 


2.5 






47909-54273 


UX Ari 


0.0 - 0.0005 


1.0 


0.00051 - 


3.5 






44928-53440 


HU Vir 


0.0-0.0001 


1.0 


0.00011 -0.0008 


1.0 


0.00081 - 


3.0 


47949-54065 


ILHya 


0.0 - 0.0005 


0.8 


0.00051 - 


3.5 






47129-54450 


YY Tri 
111 


\J.\J - yJAAAJJ 


1 n 


U.UUUJ i 










HK Lac 


0.0 - 0.0001 


0.7 


0.00011 -0.00037 


1.5 


0.00038 - 


3.2 


39751-54450 


IMPeg 


0.0 - 0.0002 


1.0 


0.00021 - 


2.0 






43734-54450 


HD 131156A 


0.0- 


1.0 










39670-52414 


HD131156B 


0.0 - 0.0005 


1.0 


0.00051 - 


5.0 






39670-52414 


HD 100180 


0.0- 


1.0 










39898-52418 


HD 201092 


0.0 - 0.00014 


1.5 


0.00015 - 0.0005 


1.0 


0.00051 - 


2.0 


39670-52518 


HD 201091 


0.0 - 0.0006 


1.0 


0.00061 - 


3.0 






39670-52518 


HD 95735 


0.0 - 0.0005 


1.0 


0.00051 - 


2.0 






39965-52419 



The values represent the different amplitudes of the cycles (see text). 

Table 5. Rotational periods from Ca II index measurements, for the full datasets and for three subintervals. 



Star 


all data 


1. part 


2. part 


3. part 




(days) 


(days) 


(days) 


(days) 


HD 131156A 


6.25 




6.21 


6.15 


HD 131156B 


11.1 




11.2 


11.05 


HD 100180 


14.3 




14.9 




HD 201092 


35.5 


30.2 


36.5 


35.53 


HD 201091 


36.1 


37.3 


36.0 


38.4 


HD 95735 


54.7 




54.5 





Values outside 3cr in different subsets indicate differential rotation. 



0.18 

X 

-^0.17 
• ^ 

^0.16 
U 

0.15 





i 

(T 


» ■■ ■ riv 


■ -tr- ^ ■ 




420 


00.0 


■ 450 


00.0 


480 


00.0 


510 


00.0 



Time-frequency distribution (STFT): 



time [J.D.] 



f [c/d] 
0.0025 



0.002 



0.0015 



0.001 



0.0005 



0.0 




42000.0 



45000.0 



48000.0 



51000.0 
time [J.D.] 



yrs 
1.10 



1.37 



1.83 



2.74 



5.48 



40000.0 



45000.0 



Time-frequency distribution (STFT); 



50000.0 

time [J.D.] 



0.001 



0.0 



40000.0 



yrs 
1.37 



1.83 



2.74 



5.48 




45000.0 



50000.0 

time [J.D.] 



